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We discuss the gravitational self-force on a particle in a black hole space-time. For a point particle, the full (bare) 
self-force diverges. It is known that the metric perturbation induced by a particle can be divided into two parts, the 
direct part (or the S part) and the tail part (or the R part), in the harmonic gauge, and the regularized self- force is 
derived from the R part which is regular and satisfies the source-free perturbed Einstein equations. In this paper, 
we consider a gauge transformation from the harmonic gauge to the Regge- Wheeler gauge in which the full metric 
perturbation can be calculated, and present a method to derive the regularized self-force for a particle in circular 
orbit around a Schwarzschild black hole in the Regge- Wheeler gauge. As a first application of this method, we then 
calculate the self-force to first post-Newtonian order. We find the correction to the total mass of the system due to 
the presence of the particle is correctly reproduced in the force at the Newtonian order. 



I. INTRODUCTION 

Thanks to recent advances in technology, an era of gravitational wave astronomy has arrived. There are already 
several large-scale laser interferometric gravitational wave detectors that are in operation in the world. Among them 
are LIGO [1], GEO-600 [2] and TAMA300 [3]. VIRGO [4] is expected to start its operation soon. The primary targets 
for these ground-based detectors are inspiralling compact binaries, and they are expected to be detected in the near 
future. 

On the other hand, there is a future space-based interferometric detector project LISA [5] that can detect gravi- 
tational waves from solar-mass compact objects orbiting supermassive black holes. There is also a future plan called 
DECIGO [6]. To extract out physical information of such binary systems from detected gravitational wave signals, 
it is essential to know the theoretical gravitational waveforms accurately. The black hole perturbation approach is 
most suited for this purpose. In this approach, one considers gravitational waves emitted by a point particle that 
represents a compact object orbiting a black hole, assuming the mass of the particle (/i) is much less than that of the 
black hole (M); /i < M. 

In the lowest order in the mass ratio {ji/MY', the orbit of the particle can be represented a geodesic on the 
background geometry of a black hole. Already in this lowest order, by combining with the assumption of adiabatic 
orbital evolution, this approach has been proved to be very powerful for evaluating general relativistic corrections to 
the gravitational waveforms, even for neutron star-neutron star (NS-NS) binaries [7]. 

In the next order, the orbit deviates from the geodesic on the black hole background because the spacetime is 
perturbed by the particle. We can interpret this deviation as the effect of the self-force on the particle itself. Since it 
is essential to take account of this deviation to predict the orbital evolution accurately, we have to derive the equation 
of motion that includes the self-force on the particle. The self-force is formally given by the tail part (or the R part 
by Detweiler and Whiting [8]) of the metric perturbation which is regular at the location of the particle. 

The gravitational self-force is, however, not easily obtainable. There are two main reasons. First, the full (bare) 
metric perturbation due to a point particle diverges at the location of the particle, hence so does the self-force. As 
mentioned above, one has to identify the R part of the metric perturbation to obtain a meaningful self-force. However, 
the R part cannot be determined locally but depends on the whole history of the particle. Therefore, one usually 
identifies the divergent part which can be evaluated locally (called the S part) to a necessary order and subtract 
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it from the full metric perturbation. This identification of the S part is sometimes called the subtraction problem. 
Second, the regularized self-force is formally defined only in the harmonic gauge because the form of the S part is 
known only in the harmonic gauge, whereas the metric perturbation of a black hole geometry can be calculated only in 
the ingoing or outgoing radiation gauge in the Kerr background, or in the Regge- Wheeler gauge in the Schwarzschild 
background. Hence, one has to find a gauge transformation to express the full metric perturbation and the divergent 
part in the same gauge. This is called the gauge problem. 

In this paper, as a first step toward a complete derivation of the gravitational self-force, we consider a particle 
orbiting a Schwarzschild black hole, and propose a method to calculate the regularized self-force by solving the 
subtraction and gauge problems simultaneously. Namely, we develop a method to regularize the self-force in the 
Regge- Wheeler gauge. The regularization is done by the "mode decomposition regularization" [9], which is effectively 
the same in the present case as the "mode-sum regularization" developed in [10 12]. 

Recently, Barack and Ori [13] proposed what they call the intermediate gauge approach to the gauge problem. 
Applying this method, the gravitational self-force for an orbit plunging straight into a Schwarzschild black hole was 
calculated by Barack and Lousto [14]. It is noted that, although their approach is philosophically quite different from 
our present approach, practically both approaches turn out to give the same result as far as the Regge- Wheeler gauge 
calculations are concerned. 

As for the case of the Kerr background, the only known gauge in which the metric perturbation can be evaluated 
is the radiation gauge formulated by Chrzanowski [15]. However, the Chrzanowski construction of the metric per- 
turbation becomes ill-defined in the neighborhood of the particle, i.e., the Einstein equations are not satisfied there 
[13]. Some progress was made by Ori [16] to obtain the correct, full metric perturbation in the Kerr background. The 
regularization parameters in the mode-sum regularization for Kerr case are calculated by Barack and Ori [17]. 

The paper is organized as follows. In Sec. II, we briefiy review the situation of the self-force problem and explain 
our strategy. In Sec. HI, we give the regularization prescription under the Regge- Wheeler gauge condition. In Sec. IV, 
we calculate the full metric perturbation and the full force in the Regge- Wheeler gauge with the Regge- Wheeler-Zerilli 
formalism. In Sec. V, we evaluate the singular, divergent part in the harmonic gauge by local analysis at the particle 
location and expand it in the Fourier-harmonic form. In Sec. VI, we calculate the S part under the Regge- Wheeler 
gauge condition by using the gauge transformation. By subtracting this S part from the full force evaluated in Sec. IV, 
we obtain the regularized gravitational self-force in Sec. VII. Finally, we summarize our calculation and discuss the 
future work in Sec. VIII. Some details of the calculations as well as discussions on the £ = and 1 modes are given 
in Appendices A F. 



II. GAUGE PROBLEM 

We consider the linearized metric perturbation 

hij.„ = Qnu - g^iy , (2.1) 

where jf^y and g^i, is the background and the perturbed metric, respectively. Here we define the force due to the 
metric perturbation as the part that gives rise to a deviation from the background geodesic: 

where z"{t) is an orbit of the particle parametrized by the background proper time (i.e., g^^idzi^ / dT){dz" / dr) = — 1). 
From the geodesic equation on we obtain 

F°'[h] = -iJ.P^(h/3^.s - ^9Dih\-s - ^Ks;l3 + ^gjsh\-0)u' , (2.3) 

where Pa^ = 6a^ + UaU^, ha/s = /la/3 — ^gaishfj,^ and = dz'^/dr. 

In the case that the perturbation is produced by a point particle, however, we face the problem that diverges at 
the location of the particle, and so does the force. Therefore, we cannot naively apply the above calculation to obtain 
the self-force of the particle. Mino, Sasaki and Tanaka [18] and Quinn and Wald [19] gave a formal answer to this 
problem by considering the metric perturbation in the harmonic gauge. According to them, the metric perturbation 
in the vicinity of the orbit can be divided into two parts: the direct part and the tail part. The direct part has 
support only on the past null cone of the field point and diverges in the limit z^{t). The tail part has 

support inside the past null cone and gives the physical self-force which is regular at the location of the particle. 
But it is almost impossible to calculate the tail part of the metric perturbation directly, because it depends on the 
global structure of the space-time as well as on the history of the particle motion. In contrast, the direct part can be 
evaluated locally in terms of geometrical quantities. Hence, instead of directly calculating the tail part, we consider 
the subtraction of the direct part from the full metric perturbation, where the latter can be calculated in principle by 
the Regge- Wheeler-Zerilli or Teukolsky formalism for black hole perturbations [20-24]. 
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From the fact that F" is a linear differential operator on h^^, (with a suitable extension of off the particle 

trajectory), we can calculate the self- force by subtracting the direct part from the full force under the harmonic gauge 
as 

lim i?^„[/i*^"'«(a;)] = lim {F^[h^^^'^{x)] - F^[h^^'^{x)]) , (2.4) 

X — >z(t) X — >z(t) 

where the superscript H stands for the harmonic gauge. When we perform this subtraction, the full metric perturbation 
and the direct part must be evaluated in the harmonic gauge because this division is meaningful only in this gauge. 

However, it is difficult to obtain the full metric perturbation directly in the harmonic gauge. In order to overcome 
this difficulty, one possibility is to perform the gauge transformation to the harmonic gauge from the gauge in which 
the full metric perturbation is obtained. In our previous paper [23] , we investigated this problem for the Schwarzschild 
case, namely, we formulated a method to perform the gauge transformation from the Rcgge- Wheeler (RW) gauge to 
the harmonic gauge. We expressed the gauge transformation equations in the Fourier-harmonic expanded form and 
derived a set of decoupled equations for the cocfhcients of each mode. Applications of this method are now under 
study. 

Recently, Detweiler and Whiting found a slight but important modification of the above division of the metric 
perturbation [8]. The new direct part, called the S part, h^'^, is constructed to be an inhomogeneous solution of the 
linearized Einstein equations (in the harmonic gauge) as 

;^full/S,H;a ^ ^j^^a ^0-^iun/S,U ^ _ ^2.5) 

The new tail part, called the R part, h^i^, is then a homogeneous solution. Since the S and R parts are both the 
solutions of the Einstein equations, wc can define the S and R parts in another gauge, which are also the solutions 
of the Einstein equations, by performing the gauge transformation of each part. Therefore, we can consider the 
subtraction procedure under some other convenient gauge by transforming the S part from the harmonic gauge to 
the desired gauge. Thus, another, perhaps more promising possibility is to formulate a method to derive the S part 
in the Regge- Wheeler or radiation gauge, where we have formalisms to evaluate the full metric perturbation, and to 
obtain the regularized self-force by subtracting the S part in this gauge. In this paper, we focus on the Schwarzschild 
case and consider the subtraction in the Regge- Wheeler gauge. 

To subtract the S part, we adopt the mode decomposition rcgularization [9]. In this method, the subtraction 
procedure (2.4) is done at each harmonic mode. The full force is obtained in the form of the Fourier-harmonic 
expansion. The Fourier (frequency) integral can be easily done in the case of circular orbits. On the other hand, the 
S part is known only in the vicinity of the particle. Hence, one has to extend it over the sphere to obtain its harmonic 
coefficients. This procedure introduces some ambiguity in the harmonic expansion of the S part. In particular, each 
harmonic mode obtained by this extension has no physical significance by itself. The physical significance is recovered 
only after we sum over all the modes. Because of this ambiguity, we have to treat the £ = and 1 modes with special 
care, as will be shown later. 



III. SELF-FORCE IN THE REGGE- WHEELER GAUGE 

The Schwarzschild metric is given in the standard Schwarzschild coordinates as 

g^^dx^dx" =-f{r)dt'^ + f{r)-^dr^ +r'^{d0'^ +sm^0d(l)^), /(r) = 1 - ^ . (3.1) 

We denote the location of the particle at its proper time r = tq as 

{^o"} = {z"{to)} = {to, ro, 00, M . (3.2) 

Formally, the gravitational self-force acting on the particle is given by the tail part in the harmonic gauge, as 
expressed in the left-hand side of Eq. (2.4). Using the notions of the S and R parts introduced by Detweiler and 
Whiting [8], it may be rewritten as 

FH(r) = \im F„[hX]ix) 

x—>z{t) 

= lha F„[h';^^'^-hl'^]{x) 

x—>z{t) 

= fim (F« (x) - F^ [hlf] {x)) , (3.3) 

x—>z{t) 

where h^'^ and h^;}^ denote the S and R parts, respectively, of the metric perturbation in the harmonic gauge. The 
S part can be calculated by the local coordinate expansion [9]. 
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Now, we consider the gauge transformation from the harmonic gauge to the RW gauge defined by 

- x^^ = -;:+e™, (3.4) 

C - hl^ = h%-'^^i.eirr (3.5) 
where is the generator of the gauge transformation. Then the self-force in the RW gauge is given by 

Ff^{T) = lim F^[h^'^'^]= lim F„ [/i^^'" - 2 V^«^^w [/i^^'H] ] (a;) 

x—*z{t) X — >z{t) 

= lim F„ - /i^.H _ 2 v^H^RW [/^fun,H _ ^s,H] ] 

x->z(t) 

= lim F„ - 2 V^H^I^W [^jfull.H] _ ^S,H ^ 2 V^H-I^W [;jS,H] ] 

= lim fF4/if"™](^)_^^[^s,H_2v^H^RW[/^s,H]](^)y (3_g) 

where we have omitted the spacetime indices of /i^j/ and V(^^j/) for notational simplicity. The full metric perturbation 
h^^u'^ can be calculated by using the Regge-Wheeler-Zerilli formalism, while the S part /i^'J^ can be obtained with 
sufficient accuracy by the local analysis near the particle location. Thus the remaining issue is if we can unambiguously 
determine the gauge transformation 

Note that the self-force (3.6) is almost identical to the expression obtained in the intermediate gauge approach [13], 
if we replace the S and R parts by the direct and tail parts, respectively. The only difference is that the S and R 
parts arc now solutions of the inhomogcneous and homogeneous Einstein equations, respectively. Hence the S part 
in the RW gauge is (at least formally) well-defined provided that the gauge transformation of the S part, Eq. (3.7) 
is unique. As will be shown later in Eqs. (6.4), this turns out to be indeed the case. Therefore one may identify the 
self-force (3.6) to be actually the one evaluated in the RW gauge [25], not in some intermediate gauge. 

IV. FULL METRIC PERTURBATION AND ITS FORCE 

In this section, we consider the full metric perturbation and its self-force in the case of a circular orbit. First, the 
metric perturbation is calculated by the Regge-Wheeler-Zerilli formalism in which a Fourier-harmonic expansion is 
used because of the symmetry of the background spacetime. Next, we derive the self-force by acting force operators 
and represent it in terms of t mode coefhcients after summing over u) and m for the Fourier-harmonic series. 

A. Regge-Wheeler-Zerilli formalism 

On the Schwarzschild background, the metric perturbation h^^ can be expanded in terms of tensor harmonics as 



Cm 



+JU{e + l)ii-l){£ + 2)Gem{t,r)fem + (V2Kem{t,r) - '^^' Gemit,r) I 



£(€ + 1), 



hoim{t,r)c\^ + — — hum{t,r)cim + ^ h2im{t,r)dem 



, (4.1) 



where a^^, aim ■ ■ ■ are the tensor harmonics introduced by Zerilli [21]. The energy- momentum tensor of a point 



particle takes the form 



J-oo ^ " dr dr 

= ^^^uV^±^^s^'Hn-n,it)), (4.2) 
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where the orbit has been expressed as 

x^'=z^'{T)={to{T),ro{T),9o{T),Mr)}, (4.3) 

with r being regarded as a function of time determined hy t = T{t). The RW gauge is defined by the conditions on 
the metric perturbation as 

The Regge- Wheeler and ZeriUi equations are obtained by plugging the metric perturbation (4.1) in the linearized 
Einstein equations and Fourier decomposing them. (Recently, the Regge- Wheeler-Zerilli formalism is improved by 
Jhingan and Tanaka [24].) 

For odd parity waves that are defined by the parity (—1)'^+^ imder the transformation (0, </>) (w — 0,(j) + tt), we 

introduce a new radial function R'l^^^r) in terms of which the two radial functions of the metric perturbation are 
expressed as 

2 

;,RW _ ^ p(odd) 

'ntmu! (^ _ ^(mu> ' 

,RW _ iAf^R(°dd)x 87rr(r - 2M) 

The new radial function R'"^^^ (r) satisfies the Regge- Wheeler equation, 

n(odd) r 2 ^r (\] p(odd) 

-^^emuj + - n[r)\R^^^ 

S.m r-2M 



[i^(£+l)(f-l)(^ + 2)]i/2 r2 

X (-^'|:[(1 - + {r- 2M)[{£ - + 2)]i/2g,„„^ , (4.6) 

where r* = r + 2Mlog(r/2M — 1), and the potential Ve is given by 



The source term Qemu> vanishes in the case of a circular orbit and 



1/2 . 0^2 

H ^-f- 5{r - ro) m de Y;^{0o, 0o) , (4.6 



where the orbit is given by 



where Cl = sjMjr^ is the orbital frequency. The orbit is assumed to be on the equatorial plane without loss of 
generality. 

For even parity waves with the parity (—1)^, we introduce a new radial function -R^^^^(r) in terms of which the four 
radial functions of the metric perturbation are expressed as 



.RW ^ A(A + \y + 3AMr + SM^ (z) r-2M d (z) 
r{r-2M)~ i{r-2Mf ~ 



Ar + 3M """" r(Ar + 3M) 
„RW _ ■ Ar^ - 3AMr - 371/^ (z) . d (z) 
tium. - *'^(^_2M)(Ar + 3M)^^™- ^^^em. 

iior^ ~ Lur{r — 2M)~ 

Ar + 3M r(Ar + 3M) 
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ttRW 

til 



Xr{r - 2M) - uj^r'^ + M{r - 3M) 



ttRW 



RW 



M(A + l)-wV3 
iujr{\r + 3M) 



rjRW 



(r-2M)(Ar + 3M) 



ttRW , fl„ 



(4.10) 



where 



X=-{i-l){i + 2), 



(4.11) 



and the source terms are given by 
87rr2(r- 2M) 



^42. + J^^m- - 167rr[i£(^ + 1)(£ - !){£ + 2)\-^I''F,^^ , 



C: 



2(.moj — — - 



Snr^ [m + 1)]-^/' 



r-2M 



r(0) 



B,mu. + -^^A^£{i+m-m+'^)r^'^Pem.. (4.12) 



r-2M 



r-2M'2 



Here the harmonic coefficients of the source terms Aimuj , ^imu '^'^'^ Bg^^ vanish in the circular case and 



-B, 



(0) 



£(£ + 1) 



-1-1/2 



Hu'^ 1 



2M\ 1 



r I r 



5{r -ro)mYl^{eo, M , 



(4.13) 



The new radial function R\J^^ (r) obeys the Zerilli equation, 



(4.14) 



where 



yf)(r) 



/ _ ^\ 2A2(A + ly + 6A^ Mr^ + 18AMV + 18M^ 
V ~) r3(Ar + 3M)2 ' 



(4.15) 



and 



(Z) 



- 2M d 
r dr 
. (r-2M)2 



(r-2M)2 / ir 



r(Ar + 3M) \r-2M 



r(Ar + 3M)2 



A(A + l)r2 + 3AMr + 6M2 - . Ar^ - 3AMr - 3M2 - 

9 L/2£mu; + « 7 ;rT7T <^Wma> 

(r — 2M) 



(4.16) 



The Zerilli equation can be transformed to the Regge- Wheeler equation by the Chandrasekhar transformation if 
desired, as shown in Appendix A. However, here, we treat the original Zerilli equation. 



B. Full Metric Perturbation 



The homogeneous solutions of the Regge- Wheeler equation are discussed in detail by Mano et al. [26] and in 
Appendix A. By constructing the retarded Green function from the homogeneous solutions with appropriate boundary 
conditions, namely, the two independent solutions with the in-going and up-going wave boundary conditions, we can 
solve the Regge- Wheeler and Zerilli equations to obtain the full metric perturbation in the RW gauge. Here, we 
consider the radial functions up to the first post-Newtonian (IPN) order. 

The radial function for the odd part of the metric perturbation metric perturbation is obtained as 

16 i n fifP mr /^^^ 



. J (2^+l)^(£+l)(£ + 2) KroJ ^^^-(^-^o) f°-<^o ' 
\2£ + l){£-l)£{£+l) I7J '^'^ 
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where O = u'^ /u*. For the even part, the radial function is obtained as 

(z) _ Snrnwu^n (( ^ -£ + A)ror 



{2£ +!){£ + 2) {£+l)uj\\ {2£ + 3)ro £{2£-l)i£ -1) 
r {£^ -2£-l) Mr 2{£'^ + £^ ~ 6£^ - 4£ - A) M 
^^n^ ^ {i-l)rl £{£-l){£ + 2)ro 

X (-^^ y;^{0o, for r < ro , 



{2£+l)£{£-l)uj\\2£-l + 1) (2£ + 3) (^ + 2) 

{£^+4£ + 2)M 2(£^ + 3£^-3£^-7£-6)M 
+ ~ i£ + 2)ro ^ (^+1) (^-1) (^ + 2)r 

X (^)V;„(0o, <^o) for r>ro, (4.18) 

The metric perturbation in the RW gauge is obtained from Eqs. (4.5) and (4.10). 



C. Full Force 

Formally, the force derived from the full metric perturbation is given by 

^f^u,Rw(-) = -f (9"- + uV) {2hit:r - C;r) ""^'^ • (4-19) 

If we decompose the above into harmonic modes, each mode becomes finite at the location of the particle though 
the sum over the modes diverges. We therefore apply the 'mode decomposition regularization' method, in which the 
force is decomposed into harmonic modes and subtract the harmonic-decomposed S part mode by mode before the 
coincidence limit x z{t) is taken. 

Since the orbit under consideration is circular, the source term contains the factor S{u — mfl), and the frequency 
integral can be trivially performed. Hence we can calculate the harmonic coefficients of the full metric perturbation 
in the time-domain. This is a great advantage of the circular orbit case, since the S part can be given only in 
the time-domain. We also note that the ^^-component of the force vanishes because of the symmetry, and F'^{z) = 
[(ro - 2M)/{r^n)]F* for a circular orbit. 

The even and odd parity parts of the full self-force are expressed in terms of the metric perturbation as 



ro (ro-3Af) (m -2M) 



^loZ = E ^^',°_3m" f'oZ,mn{ro)) do Y^Oo, <t>o) , (4.20) 



It is understood that the derivatives appearing in the above expressions are taken before the coincidence limit. It may 
be noted that there is no contribution from the components and to the even force and no contribution 

from hf^ to the odd force for a circular orbit. 

Inserting the metric perturbation under the RW gauge to the above, and performing the summation over m, we 
find 



■'^full.Rwl^ — 0' 



{£ + 1) ii'^ 1 /x^(12£3 + 25£V4£-21)M 
^ ^2 4{2£ + 3) {2£-l) 
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e ~ ~rf~2 (2£- l)(2^ + 3) ' 

^ = 0. (4.21) 

We see that the only non-vanishing component is the radial component as expected because there is no radiation 
reaction effect at IPN order. In the above, the indices (+) and (— ) denote that the coincidence limit is taken from 
outside (r > ro) of the orbit and inside (r < ro) of the orbit, respectively, and the vertical bar suffixed with £, 



full,RW 



fuU.RW 



denotes the coefficient of the £ mode in the coincidence limit. The formulas for the summation over m are shown in 
Appendix F. 

We note that the above result is valid for £ > 2. Although the £ = Q and 1 modes do not contribute to the self- force 
formally, because of our inability to know the exact form of the S part, it turns out that we do need to calculate the 
contributions from the £ = and 1 modes. These modes are treated in Appendix E. 



V. S PART OF THE METRIC PERTURBATION AND FORCE 

In this section, we calculate the S part of metric perturbation and its self- force (S-force) by using the local coordinate 
expansion. The S part of the metric perturbation in the harmonic gauge is given covariantly as 



5^a(a;, Zret)5i//3(a;, Zret)M"(Tret)u'^(Tret) 



+2/U(Tadv - Tret)fi'/:^"(a;, Z,:et)9/{x, Zret)Rja50{ZTetW {Tveth^Tret) + 0{y'^) , (5.1) 

where ^ret = -z('^ret), 7"ret IS the retarded proper time defined by the past light cone condition of the field point x, 
Tadv is the advanced proper time defined by the future light cone condition of the field point x, g^a is the parallel 
displacement bi-vector, and y is the expansion parameter of the local expansion, which may be taken to be the 
difference of the coordinates between x and zq, — x^^ — Zq . Details of the local expansion are given in [9]. The 
difference between the S part and the direct part appears in the terms of 0{y), i.e., the second term on the right-hand 
side of Eq. (5.1). In the local coordinate expansion of the S part, it is convenient to use the quantities 

e := (rg -h - 2 ro r cos 9 cos $) , 
T := i - io , R:=r-ro, 

Q:=e-^, -00. (5.2) 



A. S part of the metric perturbation 

Using the variables defined in Eqs. (5.2), it is straightforward to calculate the S part to IPN order. Here we note 
that, in general, we have to evaluate the S part up through the accuracy of 0{y), because the force is given by first 
derivatives of the metric components. The result takes the form, 

h^-^ - a T c T-R-'Q"^'' .53^ 

"'fj.v r"" / , '-m,n,p,q,r ^ ) \^-'-'J 

m,n,p,q,r 

where m, n, p, q and r are positive intergers. The explicit expressions for the components are shown in Appendix D, 
Eqs. (Dl). 



B. Tensor harmonics expansion of the S part 



In the previous subsection, we calculated the S part of the metric perturbation in the local coordinates expansion. 
In order to use them in the mode decomposition regularization, it is necessary to expand them in terms of tensor 
spherical harmonics, which involves an extension of the locally expanded S part to a quantity defined over the sphere. 
Since the only requirement is to recover the local behavior near the orbit correctly, there exists much freedom in the 
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way of extending the locally known S part to a globally defined (but only approximate) S part on the whole sphere. 
To guarantee the accuracy of /i^^ up through 0{y) in the local expansion, because the leading term diverges as 1/y, 

a spherical extension must be accurate enough to recover the behavior at 0{y^) beyond the leading order. Below, 
using one of such extensions as given in Appendix B, we derive the harmonic coefficients of the S part. 

Once we fix the method of spherical extension, it is possible in principle to calculate the harmonic coefficients of 

the extended S part exactly. However, it is neither necessary nor quite meaningful because the extension is only 
approximate. In fact, corresponding to the fact that all the terms in positive powers of y vanish in the coincidence 
limit, it is known that all the terms of 0(1/L^) or higher, where L = £ + 1/2, vanish when summed over £ [9] in the 
harmonic gauge. It should be noted, however, this result is obtained by expanding the force in the scalar spherical 
harmonics. In our present analysis, we employ the tensor spherical harmonic expansion. So, the meaning of the index 
£ is slightly different. Nevertheless, the same is found to be true. Namely, by expanding the S-part of the metric 
perturbation in the tensor spherical harmonics, the S-force in the harmonic gauge is found to have the form, 



Fi 



'M(±)| _ 



S,H 



(5.4) 



where A'' and are independent of £, and the ± denotes that the limit to tq is taken from the greater or smaller 
value of r, and 



+ 



el" 



+ ■ 



' (L2_1)(L2_4) (L2_l)(i2_4)(2,2_9) 



+ • 



(5.5) 



Then the summation of over £ (from ^ = to oo) vanishes. For convenience, let us call this the standard form. 

As wc shall sec later, the standard form of the S-forcc is foimd to persist also in the RW gauge. 

For the moment, let us assume the standard form of the S-force both in the harmonic gauge and the RW gauge. 
Then we may focus our discussion on the divergent terms. When wc calculate the S-force in the RW gauge, we 
first transform the metric perturbation from the harmonic gauge to the RW gauge, and then take appropriate linear 

combinations of their first derivatives. We then find that the harmonic coefficients h^'^, h'^im^ ^^'^ ^Um^ 
differentiated two times, and is differentiated three times, while the rest are differentiated once, to obtain the 

S-force. So, it is necessary and sufficient to perform the Taylor expansion of the harmonic coefficients up to 0{X'^) 



for/)^'" ft^"^*^'" 



and h'f^^'^, and up to 0{X^) for G^^, and the rest up to 0{X), where X = T or R. 



To the accuracy mentioned above, the harmonic coefficients of the S part are found in the form, 

'AiTmroiL'^ -2) {u'l')^ 



S,H / 



where we have defined 



r) = 


2 


r) = 


2 


r) = 


2 



£(2) (£2 _ 1) 

-2irom,{2r„ + R) {u'>')^ 



£(2) (L2 _ 1) 
lrom(72roi?L'* + 48roi?L^ + 
"6~ 



£(4) (2,2 _ 1) (^2 _ 4) 



YLi^o, 0o) , 



etc.. 



(5.6) 



£(2) = £{£+!)= i^L^- 
£(4) = ^(^+i)(^_i)(^ + 2)=(i2-i) 



(5.7) 



The explicit expressions for the coefficients are given in Appendix D, Eqs. (D2). Shown there are the coefficients in the 

case when wc approach the orbit from inside (r < tq). The results in the case of approaching from outside (r > tq) are 
obtained in the same manner. For readers' convenience, these are placed at the web page: http://www2.yukawa.kyoto- 
u-ac-jp/'misao/BEPC/. 

Now we consider the S-force in the harmonic gauge. It is noted that the t, 6 and (/(-components of the S-force vanish 
after summing over m modes. The r-component of the S-force is derived as 



S,H 



E 



2 7rAt2 



2L 



2r? 



1 M(10L^ 
- + — ^ 



11 L2 _ lOL- 17) 



4rg(L2_i) 

M (64 + 28 - 320 - 695 + 256 L -|- 442) 

16r3£(2) (L2_ i)(i2_4) 

M(156L2 _ 179) m'' 



4rg£(2) (L2-l)(L2_4)(i,2_9) 



\Yem{Oo, Ml' 
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+ 



r3£(2) (L2_i)(i2_4) 

M(2i- l)(2i2 + 2L- 1) 



The formulas for summation over m are summarized in Appendix F. For example, we have 



L 



(5i 



(5.9) 



Using these formulas, we obtain 



-^S,H 



S,H 



■ S,H 



S,H 



1 ;u2(2ro-3M) ^_ 1 /i2(4ro-7M) 

^2 4 "8 4 

M (172 - 14784 + 299) 
128rg(L2-l)(i2_4)(L2_9) 

l;u2(2ro-3M) 1/^2 (4^0 _7M) 

■ — 



+ 



^2 
0, 
0. 



(5.10) 



This is indeed of the standard form. In particular, the factor £^2) which is present in the denominators before summing 
over m turns out to be cancelled by the same factor that arises from summation over m. If it were present in the 
final result, we would not be able to conclude that the summation of over (. vanishes. We note that, apart from 
the fact that the denominator of the term takes the standard form, the numerical coefficients appearing in the 
numerator should not be taken rigorously. This is because our calciilation is accurate only to 0{if') of the S-force, 
while the numerical coefficients depend on the 0{y) behavior of it (an example is shown in Appendix C). It is also 
noted that the 0(l/X)-terms are absent in the S-force, implying the absence of logarithmic divergence. 

It is important to note that i in the above runs from to oo. Although there are some tensor harmonics that do 
not exist for 1 = and/or C. = 1, we may regard that the corresponding harmonic coefficients contribute to the 
and £>^ terms of the S-force individually, with + D'^ = 0. That is, we set the contributions to A^^ to zero and 
adjust the term in such a way that D'^ = —B^ for these special coefficients while keeping the standard form for 



VI. S PART IN THE REGGE- WHEELER GAUGE 

Now, we transform the S part of the metric perturbation from the harmonic gauge to the RW gauge. The gauge 
transformation functions are given in the tensor-harmonic expansion form as 

^(odd) 



em 



^(even) ^ ^|Mo«^™(t, r)r,„(0, 0), M,^ 

Mll^-^^'^it, r)deYemiO, </>) , M^,^^-^^ (t, r)9^F,„(0, 0) \ . (6.1) 



There are one degree of gauge freedom for the odd part and three for the even part. To satisfy the RW gauge 
condition (4.4), we obtain the equations for the gauge functions that are found to be rather simple: 

hlfjt,r) = -2.AL™(t,r), 
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h, 



(e)S,H 
Mm 



fS.H^RW 
2lm 



^ ,^S,H^RW 



(t, r). 



(6.2) 



We therefore find 



(6.3) 



^■™(i,r) = -^GL«(i,r) 



2^m 



Mo^i™(t,r) = -4r(t,r)-a,M 



Ml (t, r) 



(e)S,H 



S.H^RW 
'■21m 



it. r) ., 



it, r) - a. 



M 



S.H^KW 
2im 



it, r) 



(6.4) 



Wc note that it is not necessary to calculate any integration with respect to t or r. It is also noted that the gauge 
functions are determined uniquely. This is because the RW gauge is a gauge in which there is no residual gauge 
freedom (for £ > 2). 

Then the S part of the metric perturbation in the RW gauge is expressed in terms of those in the harmonic gauge 
as follows. The odd parity components are found as 



Kern it, r) 



Kimit^n+dtly^'^ it,r), 

^tm it' n 



(6.5) 



and the even parity components are found as 



o-S.RW/, X _ 
-"O^m V^' ^> — 



K, 



S,RW 
Im 



2M 



Um ^) 



HfiTit, r) = HlZit, r) 



S,H^RW 



it, r) + drM^l"it, r) 



2M 



-m: 



S.H^RW 



M^l^ it, r) + 



r(r -2M) 
r(r - 2 M) "™ 



oem 



it, r) 



it, r) 



it, r) = K^e^it, r) + 



2(r-2M) ^^g,H 



+RW 



Um 



it, r) , 



(6.6) 



where the gauge functions are given by Eqs. (6.3) and (6.4). 



A. Gauge transformation and the S psirt in the RW gauge 



Inserting the results obtained in the previous section to Eqs. (6.3) and (6.4), we obtain the gauge functions that 

transform the S part from the harmonic gauge to the RW gauge. They are shown in Appendix D, Eqs. (D3). It may 
be noted that the gauge functions do not contribute to the metric at the Newtonian order. In other words, both the 
RW gauge and the harmonic gauge reduce to the same (Newtonian) gauge in the Newtonian limit. 
The S part of the metric perturbation in the RW gauge is now found in the form. 



h: 



S,RW 

oem 



it, r) = 



h'^nt' r) 



TT jJL 



AiTmroiL^ -2) iu'^'f 



+ ■ 



£(2) (L2 _ 1) 
imro (-60roi^ + 174roi^ + • 



9,y;„(^o, </'o), 
)(u*)2i 



£(4) (2,2 _ 1) (i2 _ 4) 



Yemi^O, ^o) , 



etc. 



(6.7) 



The explicit expressions are given in Appendix D, Eqs. (D4). 
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B. S force 

Next we calculate the S part of the self-force. Of course, it diverges in the coincidence limit. However, as we noted 
several times, in the mode decomposition regularization in which the regularization is done for each harmonic mode, 
the harmonic coefficients of the S part are finite. 

The calculation is straightforward. Expanding the formula for the self- force (2.3) in terms of the tensor harmonics, 
we obtain 



pf RW 
(even) 



= E 



F, 



rRW 
(even) 



ptRW 



■^(odd) 



em 



2(ro-3M)2(ro-2M) 



-ro M {dt H^Zih, ro)) + 2 {d, H^Zih, ro)) 



E 



+imrlnH^Z{to, ro)-6imMroflH^Z{to, ro) + 8imM^nH^Z{to, ro) 
-imroMQ Kf^{to, ro) + 2imM^n Kf^{to, ro) + 5 {d^ Kf^{to, ro)) 

-2roM{dtKf^{to, ro))) F^„(^o, "Ao) , 

2 rl {dt HfZito, ro)) - 2MH^Z(.to, ro) 



I2{>u-2M) 



em 2r§(ro-3M) 

+2imrlft HfZito, ro) - 2 M Kf^{to, ro) - r^ {dr H^Zito, ro)) 
+2 ro M {dr H^Zito, ro)) - ro M (9, Kf^{to, ro))) F^„(^o, M , 



E - 



em 



(ro-3M)2(ro-2M) 



nmMhX^{to, ro)-i{ro-2M) {dth^Z^i^o, ro))) 



xdgYim{do, M , 

{{drhTm{to,ro))- {dthfrm{to,ro))-inmhfZ{to,ro)) 



em 



ro-3M 
xd0Yem{Oo, <f>o)- 



(6.8) 



Substituting the S part of the metric components in the RW gauge as shown in Eqs. (6.7), given explicitly in Eqs. (D4), 
into the above, we find that the t, and ^-components of the S-force vanish after summing over m. The r-component 
of the S-force inside the particle trajectory is derived as 



S,RW 



= E 



2 7rA«^ [/2L- 1 _^ M(10L3 + 11L2 - lOL- 17) 



2rl 



4r'oiL'-l) 



M (64 + 28 - 320 - 695 + 256 L + 442) 
16 r^ £(2) (L2-l)(L2_4) 

M (156^2 - 179) m-* 



+ 



Summing the above over m, we obtain 



4r3£(2) (L2-l)(L2_4)(i2_9) 
13Mto2 



\Yem{Oo, <Ao)r 



rg£(2) (i2-l)(L2_4) 

M(2L-1)(2L2 + 2L-1) 



rg£(2) (L2-1) 



\deYe^{0o, 0o)|' 



(6.9) 



p* I 

^S,RW 



-^S,RW 



S,RW 



0, 

l/x2(2ro-3M) 1/^2 (4^0 _7M) 



+ 



iLt2 M (172 - 14784 L'^ + 299) 
128rg(L2-i)(i2_4)(i2_9) 



IfU^ro^SM) lMl(4ro^7M) 1 

^2 4 ^~8 ^ +^(;^)^ 

0, 
0. 



(6.10) 
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We now sec that the S-force in the RW gauge also has the standard form as in the case of the harmonic gauge and 
there is no 0{1/L) term. Note that, again with the same reason as we explained at the end of the previous section, 
the final formula above should be regarded as valid for all £ from to oo. 



VII. REGULARIZED GRAVITATIONAL SELF-FORCE 



In the previous two sections, we have calculated the full and S parts of the self-force in the RW gauge. Now we are 
ready to evaluate the regularized self-force. But there is one more issue to be discussed, namely, the treatment of the 
£ = and 1 modes. 

The full metric perturbation and its self-force arc derived by the Regge- Wheeler- Zcrilli formalism. This means they 
contain only the harmonic modes with £ > 2. If we could know the exact S part, then the knowledge of the modes 
i>2 would be sufficient to derive the regular, R part of the self-force, because the R part of the metric perturbation is 
known to satisfy the homogeneous Einstein equations [8] , and because there are no non- trivial homogeneous solutions 
in the i = and 1 modes. To be more precise, apart of the gauge modes that are always present, the £ = homogeneous 
solution corresponds to a shift of the black hole mass and the £ = 1 odd parity to adding a small angular momentum 
to the black hole, both of which should be put to zero in the absence of an orbiting particle. As for the £ ~ 1 even 
mode, it is a pure gauge that corresponds to a dipolar shift of the coordinates. In other words, apart from possible 
gauge mode contributions, the £ — and 1 modes of the full force should be exactly cancelled by those of the S 
part. In reality, however, what we have in hand is only an approximate S part. In particular, its individual harmonic 
coefficients do not have physical meaning. Let us denote the harmonic coefficients of the approximate S-force by 

F^''^^, while the exact S-force and the full force by Ff and J^/"", respectively. Then the R-force F^ may be expressed 
as 



F^ = j2 (Fr - f!) 

e>2 

e>o e>o 

= E {Fr - ^■^) + E {Fr - f!'^^) , (7.1) 

e>2 (=0,1 

where = Ff — F^''^^, and the last line follows from the fact that F^'^^ are assumed to be obtained from a 
sufficiently accurate spherical extension of the local behavior of the S part to guarantee X]£>o ~ ^- Thus, it is 
necessary to evaluate the £ = and 1 modes of the full force to evaluate the self-force correctly. 

First, we consider the contributions oi £ > 2 to the self-force. As noted before, for the IPN calculation, the only 
r-component of the full and S part of the self-force is non-zero. The £ mode coefficients corresponding to the first 
term in the last line of Eq. (7.1) are derived as 

I rpr I Tpr \ 

^RW\i — -f^full.RWl^ ~ -f^S.RWl^ 



8(2^- l)(2f -1-3) rg 



Summing over £ modes, we obtain 



'^0 



F^y,{£>2) = -^^. (7.3) 



Next, we consider the £ = Q and 1 modes. Detailed analyses are given in Appendix E. It is noted that the £ = Q 
and £ = 1 odd modes, which describe the perturbation in the total mass and angular momentum, respectively, of 
the system due to the presence of the particle, are determinable in the harmonic gauge, with the retarded boundary 
condition. On the other hand, we were unable to solve for the £ = 1 even mode in the harmonic gauge. Since it is 
locally a gauge mode describing a shift of the center of mass coordinates, this gives rise to an ambiguity in the final 
result of the self-force. Nevertheless, we were able to resolve this ambiguity at Newtonian order, and hence to obtain 
an unambiguous interpretation of the resulting self-force. 

The correction to the regularized self-force that arises from the £ = Q and 1 modes, corresponding to the second 
term in the last line of Eq. (7.1), is found as 

^0 
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Finally, adding Eqs. (7.3) and (7.4), we obtain the regularized gravitational self-force to the IPN order as 

-^RW - — 2 3 • V'-oj 

'0 'O 

Since there will be no effect of the gravitational radiation at the IPN order, i.e., the t— and 0— components are zero, 
the above force describes the correction to the radius of the orbit that deviates from the geodesic on the unperturbed 

background. It is noted that the first term proportional to /j,^ is just the correction to the total mass of the system 
at the Newtonian order, where tq is interpreted as the distance from the center of mass of the system to the particle. 



VIII. CONCLUSION 



In this paper, we proposed a new method to derive the regularized gravitational self-force on a point particle in 
circular orbit around the Schwarzschild black hole, and, as a demonstration, we derived the regularized self-force 
analytically to the first post-Newtonian (IPN) order. 

The regularization of the gravitational self-force may be divided to the two problems, the subtraction problem and 
the gauge problem. To regularize and subtract the divergent part, we employed the 'mode decomposition regulariza- 
tion', in which everything is expanded in the spherical harmonics and the regularization is performed at each i mode. 
As for the gauge problem, utilizing the recent discovery by Detweiler and Whiting that the regularized force may be 
derived from the II part of the metric perturbation that satisfies the source-free Einstein equations, we considered the 
regularized force in the Regge- Wheeler gauge. 

In the present paper, actual calculations were done only for circular orbit and to the IPN order. However, there 
remains a problem for the even parity £ = 1 mode. In this metric perturbation approach, there inevitably remains 
ambiguity of the gauge in the resulting self-force. To circumvent this difficulty, the only way seems to be to regularize 
at the level of the Weyl scalar ip4 or the Hertz potential which are free from the i = and 1 modes As another 
problem, to make our method applicable to general cases, it is therefore necessary to extend to general orbits and 
to higher PN orders. Some progress in this direction based on analytical methods is under way [27]. It will also be 
necessary to incorporate numerical techniques if we are to treat completely general orbits. Some development is done 
by Fujita et al. [28]. 

Our final goal is to derive the self-force on the Kerr background. Recently, Mino [29] has proposed a new approach 
to the radiation reaction problem by using the radiative Green function. In his method, assuming the validity of 
the adiabatic approximation, the radiation reaction to the conserved quantities including the Carter constant can be 
calculated from the radiative Green function, which is free from any singular, divergent behavior. This is a great 
computational advantage. However, this method cannot treat the self-force for completely general orbit because of 
the assumption of adiabaticity. It is therefore still necessary to derive the self-force in the general case. One possibility 
is to consider the regularization of the Weyl scalar ^'4 and construct the R part of the metric perturbation in the 
radiation gauge by using the Chrzanowski method. Investigations in this direction is also in progress [30]. 
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APPENDIX A: MANO ET AL. ANALYSIS 

In this Appendix, we summarize the analysis of Mano et al. [26] which we use in order to derive the full metric 
perturbation for ^ > 2 modes. 



1. Homogeneous solutions 



We investigate the analytic expression of the Regge- Wheeler functions, and generate these functions in an explicit 
manner up to 0(v^) corrections relative to the leading order in the slow-motion expansion, i.e., first post-Newtonian 
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order. (More detail analysis is given in [31].) Here u is a characteristic velocity of the particle. The Regge- Wheeler 
equation is 



d_ 

dr 



2M 
r 



dr 



2M 
r 



(co'-Veir)) 



2M 
r 



<^(even/odd) 



(r). 



(Al) 



The source term Sj^^^^ is expressed in terms of the source terms of the Zerilli equations [21] as 



"((even) 



A(A + 1) + 



9M' 



r2(Ar + 3M) J 



r J dr 



and the Zerilli function, -R^^^^ is derived from -R^^™^ as 



d(z) ^ I 

^em^ (A2(A + l)2 + 9a;2M2) 



A(A+1) + 



'ir-2M) \ ^(even) 



9M 



r2(Ar + 3M) J 



'imuj 



+ 3M 1 - 



2M \ _d ^(even) 

r J dr 



(A2) 



(A3) 



So, we may focus on the Regge- Wheeler function. The Regge- Wheeler equation is rewritten as 



z — e 



2e 

1 + + 



i{e + 1) 



+ ■ 



3e 



z — e {z — e)'^ z{z — e) z'^{z — e) 



X{z) = (l-^y's{z)iA4) 



Here z = Lor and e = 2Muj, and wc use the symbol X{z) for Ri^^^°'^^\r), S{z) for Sf°^^"^^°'^'^\r). In the post- 
Newtonian expansion, both z and e arc assumed to be small, while only e is considered to be small in the post- 
Minkowskian expansion. We note that z ~ 0{v) and e ~ 0{v^) in the post-Newtonian expansion. 

First, we consider a homogeneous Regge- Wheeler function in the form of a series of the Coulomb wave functions, 
Xc". (See (3.4) and (3.6) in Ref. [26].) 



.^T{n + iy-l- ie)T{n + 



Fn+.{z) = e-"{2z 



.n+u T{n + iy + l + ie) 
T{2n + 2i^ + 2) 



r{n + iy + l + ie)T{n + i' + 3 + ie) 

iFi{n + iy + l + ie; 2n + 2v + 2; 2iz) , 



(A5) 



where i^i is the confluent hypergeometric function, and the expansion coefficients are determined by the three- 
term recurrence relation, (See (2.5) and below in Ref. [26].) 

{n + u — 1 + ie){n + ly — 1 — ie){n + p + 1 — ie) 



a. 



-le- 



(n + z/+l)(2n + 2z/ + 3) 



^ _ . {n + i/ + 2 + ie)(n + i/ + 2 - ie){n + v + ie) 
^" {n + v){2n + 2v-l) ' 



(A6) 



and V, which is called the renormalized angular momentum, is determined by requiring the convergence of the series 
expansion in Xc'^ . Replacing v by —v — 1, one obtains the other independent solution Xc" ■ It is important to note 
that the renormalized angular momentum in the post-Minkowskian expansion becomes 



(A7) 



Hence v = £ to IPN order. 

The post-Minkowskian expansion of the coefficients an" is also discussed in Ref. [26]. With the normalization 
aa" = 1, they are found for £ > 2, 



a-t 

V 

an 
an'' 



0(el"l) 
0(6^+i), 
0(e^+2), 
0(e^+^), 



{n>-i + 2), 



{-e-2>n> -2i), 
{-2£ -l>n). 



(A8) 
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The post-Minkowskian expansion of the coefficients a„ " ^ can be obtained by using the symmetry, 

= a-n-"-^ ■ (A9) 

(See (2.13) in Rcf.[26].) 

The leading terms in the Regge- Wheeler functions in the slow-motion expansion become 

Xc" ~ 0(^^+^e°), 
Xc""-^ ~ 0{z-^e°). (AlO) 

Then, for instance, if we consider IPN order, it is sufficient to take account of the Oq and a''_i terms in Xq" and 

Xc-"-^ {l>2). 

In Rcf. [26], the homogeneous Rcggc- Wheeler functions with the in-going and up-going boundary conditions are 

derived in the form of linear combinations of Xc'^ and Xc^'''^^ ■ The in-going boundary condition is that waves are 
purely in-going at the black hole horizon, and the up-going boundary condition is that waves are purely out-going at 
the infinity. 



Xin" = Ki.Xc" + K_^_iXc~ 



7r2-'"e 



, — — 1 



r{u + 1 + ie)r{i^ - 1 + ze)r(i/ + 3 + ze) sin7r(z/ + ie) 



n=0 



T{n + v -l + ie)T{n + 2y+l) ^ 



E 



g2j7i-i/ I si£iL(idi£) I sin7r(iy - ie) 

sin7r(i/— ze) 

The leading order of and K-^-i for £ > 2 becomes 



n!r(n + v + ^~i€) 

T{n + v — 1 ~ ie)T{n + 1^+1 — 16) ^ 
{-n)\r{n + 1^+1 + ie)T{n + 1^ + 3 + ie)r{n + 2i^ + 2) 

sin7r(. + ze)^^,_.^,,,^_._i- 



Then we find 



K^Xc" 

Therefore, we may replace Xi^" by Xc" to 3PN order. As for X-ap" , we find 

sin Tt{v+ie.) u 

-ie^'^'Xc'"'^ V y V / 

Thus, we may replace X^p"^ by Xc~''~^ to 2PN order. 

For convenience, we define the homogeneous solutions Xc" and X^'^~^ normalized as 



r(z/ + 3 + ie)r(2z/ + 2) 
(1 + Oiv)) , 



- r(-z.-2-fe)r(-^.-ie)^^ 
= z-'{l + 0{v)) . 



These are expanded to 0{v^) as 



Xc'-iz) = z{2zr(l-l 



1 {e-2){£ + 2)e 



2 2, + 21 2 



Iz 



Xc-"-\z) = z{2z)-''-^ 1 + 



1 z' 



+ 



1 (£- l)(£ + 3)e 



22^-1 2 {i + l)z 



+ 0(t;') 



(All) 
(A12) 



(A13) 
(A14) 

(A15) 



(A16) 
(A17) 

(A18) 
(A19) 
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where V = (. + O(f^). 

To summarize, for the in-going homogeneous sohition normahzed as 

Xin'^(z) = Xc'' + P.Xc-'-^ (A20) 
all the coefficients (3^ can be set to zero up through 3PN order, while, for the up-going solution normalized as 

X„p-(^) = Xc-'-^+l.Xc'', (A21) 
we may put 7^ = up through 2PN order. 

2. Retarded Green function 

Using the Regge- Wheeler functions X-m" and X^^" , the retarded Green function is constructed as 

Gret''(^, Z') = W{X,^'^{XX^^'^{Z')) (^in''(^)^up^(^')^(^' " ^) + X^^^ {z)X,^'^ {z')e{z - z')) , (A22) 

where W{Xi,X2) is the Wronskian, 

WiX,iz'),X2{z')) = - (1 - J) (x,{z')-^X2{z') - X2iz')-^Xi{z')^ = const. . (A23) 



This Green function satisfies 
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= -(1-^) ^5[z-z'). (A24) 
Then the Regge- Wheeler function with the source term S^^^'^^°'^'^\r) is given by 

^(ewodd)^^^ ^ rf,.G,et''(r, r')i (1 - ^) stZ'°'^\r') . (A25) 

J2M <^ V ^ / 

Here wc arc only interested in the Green function accurate to IPN order. A numerical method to construct an 
accurate Green function based on this Mano-Suzuki-Takasugi method is discussed in Ref. [28]. 

APPENDIX B: SPHERICAL EXTENSION OF S PART 

In this Appendix, we consider the tensor harmonic decomposition of the S part. First, we give the decomposition 
of e" where 

e= {r^ + 7-2 - 2ror cos U ■ Ua) , (Bl) 
and rio is taken to be on the equatorial plane, {n/2, 0o). Extending ft over the whole sphere, we have 



= E WTi^'"^^'' ro)YU^)Y;^{Slo) , (B2) 



where the detail derivation as well as the coefficients E^^'' are given in Appendix D of [9] . 
In terms of the coefficients Ei^^ , the formulas needed to decompose the S part are derived as 
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Note that these formulas are valid only in the sense of the spherical extension given by Eq. (B2). 



APPENDIX C: 0(y^)-CORRECTION 

In this Appendix, as an example to clarify how the standard form is recovered and why it is necessary to include 
the t = Q,l modes even if some of the tensor harmonics are identically zero for these modes, we consider the S part of 

the metric components hte and htrf, and analyze the contribution of their 0{y'^) terms to the self-force in the harmonic 
gauge. 

These metric components give rise to the coefficient h^^}^ of the vector harmonic proportional to {dgYfrn, d^Y^rn)- 
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Note that this vanishes identically for £ = 0. Since the contribution of the 0{y'^) terms to the self-force is zero, the 

sum of /iQ°j„ over £ should vanish. We show that it indeed has the standard form for general £. However, to guarantee 

that the sum over i is zero, it is necessary to include the contribution from £ = as well. This implies B^^ + D'^ = 

for ^ = as discussed at the end of Section V. 

The local expansion of the S part of the metric components hte and /it^ takes the form 
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where we have retained only terms that may contribute to the self- force, and the superscript means 0{y^). The 
tensor harmonic coefficients h^j^ are given by 
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where £ 0. For the 0{y'^) terms of the form (CI), we have 



The force is given by 
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Here, since the terms of interest arc already of 0{y^), we may use the leading order formulas for the spherical 
extension of the local coordinate expansion [9]. We have 
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where n > 1, r> = max{r, Tq}, r< = min{r, ro}, L = £ + 1/2 and k„ is a constant independent of L. Therefore, 
Eq. (C3) is evaluated as 
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Using Eq. (C4), and retaining only the terms that will remain after summing over m, we have 
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Thus, the O(y^) terms contribute to the term in the form of the standard form, and the sum over £ vanishes 

provided we include the £ = term in the summation. Since the 0{'tp') terms do not contribute to the force anyway, 
it then follows that we may adjust the numerators of the D'^ term so as to give Dq = — iJ^. 



APPENDIX D: CALCULATION OF THE S PART 



In this Appendix, wc show the S part of tlic metric perturbation and its gauge transformaton. The S part of the 
metric perturbation under the harmonic gauge are obtained in the local coordinate expansion as 
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The harmonic coefBcients of the above S part are calculated as 
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The gauge transformation from the harmonic gauge to the RW gauge is given by 
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And then, the coefficients of the S part under the RW gauge are calculated as 
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APPENDIX E: £ = AND 1 MODES 



In this Appendix, we derive the contributions to the self-force in the RW gauge from the £ = and 1 modes. As 
discussed at the beginning of Sec. VII, and described in Eq. (7.1), although there is no physical contribution from 
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the £ = and £ = 1 modes to the self-force in the rigorous sense, since we can calculate the S part only locally in the 
vicinity of the particle, its spherical extension inevitably contaminates each £ mode with other £ modes. Therefore, 
in particular, we have to take account of the corrections from the £ = and £ = 1 modes to the self-forcx;. 

For the £ = and £ = 1 odd modes, the RW gauge condition is automatically satisfied, since = h['^ = G = for 
£ = and /i2 = for ^ = 1 odd modes. An appropriate choice of gauge is then to consider the perturbation under the 
retarded causal boundary condition in the harmonic gauge. In fact, if we recall the gauge transformation equations 
from the harmonic gauge to the RW gauge given by Eq. (6.2), we see that all the gauge transformation generators 
for £ = and £ = I odd modes vanish. Thus, no gauge transformation is needed for the S part of these modes, and 
our task is to find the exact solutions in the harmonic gauge with the retarded boundary condition and perform the 
subtraction of the S part under the harmonic gauge. 

For the £ = 1 even mode, the RW gauge condition is non-trivial and there is a gauge degree of freedom in the 
RW gauge, reflecting the fact that it is a pure gauge mode describes a shift of the center of mass coordinates in the 
source-free case. On the other hand, the gauge transformation of this mode from the harmonic gauge to the RW gauge 
is uniquely determined. Thus, to determine the self-force unambiguously in the RW gauge, one first has to solve the 
perturbation equations in the harmonic gauge exactly (under the retarded boundary condition), transform the result 
to the RW gauge, and perform the subtraction of the S part. However, unfortunately, we were unable to solve for 
the £ = 1 even mode in the harmonic gauge due to a complicated structure of the perturbation equations (i.e., in the 
form of coupled hyperbolic equations). Thus, there remains a gauge ambiguity in the final result. Nevertheless, in the 
Newton limit when the coordinates can be defined globally, we can resolve the gauge ambiguity and give a definite 
meaning to the resulting self-force. 

To summarize, the regularized self-force in the RW gauge is expressed as 

^R,RW ^ ^ l^^fall.RW _ j,S,RW.Ap^ ^ ^ ^ (El) 

where 
with 

rpRW _ zpfull.H ,;,S.H,Ap 

<'-f^^=l(odd) - -^^=l(odd) ~ ^t. 



l(odd) ~ -^^=l(odd) ^ £=l(odd) 
l(even) " -'^^=l(even) -'^^=l(even) 

where there remains a gauge ambiguity in ^F^^^^^^^^^ 



rpRW _ pfull.RW pS,RW,Ap /t-i^n 



1. i — mode 



First, we consider the £ = Q mode of the full metric perturbation. It is noted that the £ = Q mode consists of only 
the even parity part and all the derivatives of Yoo vanish. As noted above, this mode satisfies the RW gauge condition 

(e) (e) 

h/Q = h\ = G = Q automatically. So, the appropriate choice of gauge is the harmonic gauge under the retarded 
boundary condition. To find the exact solution in this gauge, we consider a gauge transformation of the exact solution 
foimd by Zerilli. 

This mode represents the perturbation in the total mass of the system and was analyzed by Zerilli. For the £ = Q 
mode, there are two gauge degrees of freedom. The choice made by Zerilli is 



Js:f"M'Z(t,r) = 0. 



(E4) 



which we call the Zerilli (Z) gauge and denote the quantities in it by the superscript Z. In the case of a circular orbit, 
the £ = Q mode metric perturbation is solved to be 



rTfull,Z/, N 
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e(r-ro). 



ro - 2M 
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(E5) 



Here we imposed the boundary condition that the black hole mass is unperturbed and the perturbation satisfies the 
asymptotic flatness. Note that the Zerilli gauge is singular in the sense that the metric has a discontinuity at r = ro- 
The constant a is given by 



2(47r)i/V«* 



1 - 



2 M 
ro 



(E6) 



27 



Note that the £ = mode is independent of time. So we may write iJ2""'^(i, r) = H^^^''^{r). 

Now we consider the gauge transformation from the above Z gauge to the harmonic gauge. The equations for the 
gauge transformation are formally written as 



Detailed discussions on the gauge transformation to the harmonic gauge are given in [23]. 
We set the gauge transformation generator as 

U^^"} = {Mo(r)yoo(^, <^),Mi(r)Foo(^, <^),0,0} . 
In the circular case, the i = Q mode of Eq. (E7) is explicitly written down as 

" (P 2d_' 

dr"^ r dr 
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where 
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Since Mq is independent of the source, we set it to zero in accordance with the retarded boundary condition. Thus 
we focus on the equation for Mi. 

We employ the Green function method to solve Eq. (ElO). Two independent homogeneous solutions are easily 
obtained as 



^(homo,l) _ 

1 r(r-2M)' 



(homo, 2) 



r-2M 
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Using the above homogeneous solutions, we construct the other two independent solutions M™ and M°"* which are 
regular at the event horizon and infinity, respectively. We find 



+ 2Mr + 4M2 



r(r - 2M) 
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Then the Green function is derived as 

1 



and Ml is given by 



G{r,r') = — [Mj"(r)M°"^(r')e(r' - r) + Af°"*(r)Mj"(r')e(r - r')] 
W = (r-2M)2[Mi"(r)a^Mr*(r)-M°^*(r)a^Mi"(r)] =-3, 



Mi(r)= / G{r,r')r'{r' -2M)S{r')dr' . 
J2M 
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Although the integral can be performed without any approximation, we only show the result to IPN order, 

5a r a(13ro + 6r)M" 
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The metric perturbation transforms under the above gauge transformation as 
H„«(r) = /fJ(r) + i^M,(r), 



(E17) 



Note that we have Hi =0 because Mq = 0. (It may be noted that Hi does not contribute to the force for a circular 
orbit even if it is non-zero.) Then the metric perturbation in the harmonic gauge is found as 
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and the full force is calculated as 



^mi,h(^ = 0) = 
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2r3 



©(r-ro). 
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Next, wc consider the S part of the metric perturbation. Its harmonic coefficients are given in Eqs. (5.6). Only the 
harmonic coefficients Hq, H2 and K remain for the £ = mode. To IPN order, we have 
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The S-force in the harmonic gauge is calculated as 

^s.h(^ = 0) = [^^1 0(^0 -r) + 
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Prom the above results, we obtain the contribution of the £ = force as 

6F^^{£ = 0) = 5F^{e = 0) 

= ,,h(^ = 0) - Fl^{£ = 0) 
23m^M 



(E22) 



2. £. — 1 odd parity mode 



The £ = 1 odd mode represents the angular momentum perturbation added to the system. It also satisfies the odd 
parity RW gauge condition /12 = automatically. Therefore, as in the £ = case, we look for the exact solution in 
the harmonic gauge with the retarded boundary condition. 



29 



The full metric perturbation consists of the two components, /iq"'' and ft,^"''. These were also solved by Zerilli. There 
is one gauge degree of freedom, and we may put hi = 0. The appropriate boundary condition is that the black hole 
angular momentum is unperturbed and the perturbation is well-behaved at infinity. Then we find 



where b is given by 
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Note that only the m = mode is non-zero, and it is time- independent. 
Next, we consider the gauge transformation to the harmonic gauge. We set 



The equation for A^~*^ becomes 
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This is a source-free hyperbolic equation. So, with the retarded boundary condition, we find A^j~*^ = 0, that is, 
the Zerilli gauge is equivalent to the harmonic gauge with the retarded boundary condition. The full force is then 
calculated as 



-^r{odd) / 



The harmonic coefficients of the S part are given as 
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and the S-force is obtained as 

Subtracting the S part from the full force, we find 
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= 0. 



(E30) 



Thus, our spherical extension turns out to be accurate enough to reproduce the correct £ = 1 odd mode up to IPN 
order. 



3. £ = 1 even parity mode 



The £ = 1 even mode represents essentially a gauge mode that describes a shift of the center of momentum of the 
system. The coefficient G is absent from the beginning, while there is no loss in the gauge freedom. Hence there 
remains one degree of gauge freedom in the RW gauge. As mentioned at the beginning of this Appendix, to fix 
the gauge completely, it is necessary to solve the perturbation equations in the harmonic gauge with the retarded 
boundary condition, and to perform the gauge transformation to the RW gauge. However, because the perturbation 
equations become complicated, coupled hyperbolic equations in the harmonic gauge, we were unable to solve for this 
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mode. Here, we therefore give up fixing the gauge unambiguously, but solve the perturbation equations in the RW 
gauge, imposing a gauge condition by hand. 

To look for an exact solution in the RW gauge, following Zerilli, we choose K = in addition to h^Q ^ = h^^ = 0. 
Let us also call it the Zerilli gauge. The harmonic coefficients for the full metric perturbation in the Zerilli gauge are 
given by 
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where 
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(E32) 



and we have imposed the boundary condition that the perturbation is regular at horizon. It may be noted that 
although the £ = 1 even mode is locally a pure gauge, it is not so in the global sense because of the regularity at the 
horizon. Note that the m = components vanish because the orbit is on the equatorial plane. It is also noted that 

-(M/r3)/„ for 
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The coefficient ^^q"'^^ in the above behaves as ~ r at infinity. Without violating the RW gauge condition, it is 
possible to remove this singular behavior. Namely, we consider a gauge transformation, 
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As a solution of the above gauge equations that makes the metric perturbation regular at infinity, we choose 
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By the above gauge transformation, the r-component of the force changes by 
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So, to IPN order, we find 
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Thus, the full force in this RW gauge is given by 
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It may be noted that, at Newtonian order, the r coordinate of the Zerilli gauge, in which the metric inside the orbit 
is unperturbed, corresponds to placing the black hole at r = 0, while the gauge transformation that regularizes the 
asymptotic behavior at infinity makes r the radial coordinate measured in the center of mass coordinate system. In 
other words, ro in the Zerilli gauge gives the relative distance between the black hole and the particle, while ro after 
the transformation gives the distance from the center of mass to the particle. This explains the Newtonian part of 
the change in the force, 2>\j? jr^. In this sense, the gauge freedom is under control at Newtonian order. 

Now we turn to the S part. The harmonic coefficients in the harmonic gauge are given by 
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Wc transform the above to the RW gauge, as discussed in Section V. Since G is absent from the beginning, Eqs.(6.4) 
that give the gauge transformation from the harmonic gauge to the RW gauge are simphfied as 
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The resulting harmonic coefficients in the RW gauge are expressed as those given in Eqs. (6.6) except for the gauge 
functions Mq and Mi that are now given by the above equations. From these, we find 



H. 



S,RW 
11m 



{t, r) 



^^21™ r) 



Kfnt, r) 



2ro + 2R 3MR_2iTmu'f 



ro 



+ ^ {-Urom^ +9 R- 31 Rm^ + ^ro) {u'^'f 
2ro-4i? 3MR 2iTmu^ 



e(ro - r) 



''0 



2 

+ 9 



63 



27 



— R-Urom^ + 17 Rm^ + — ro] (uf) 



e(r - ro) 



TM 8 



'0 



9 



8 /9 



- {2ro + -R] imu'^-- - + 2m^ {u'^yT 



e(ro - r) 



TM 8 



„3 



9 



- [Aro- -R] imu't' - - --4mM (u'f')^ T 



9 



e(r - ro) 



4 



■ TT/U 



2ro + 2E 3MR_2iTmuf 
o H ^ 



'0 



'0 



ro 



+ 



(9i?- 2rom^ + 15 Rm^ + ^ ro) {u'^'f 
2ro-4:R 3MR 2iTmu'>' 



ro r-g ro 

ro - 2 ro - ^ i? + 15 i?m2) {u^f 



e(ro - r) 



e(r - ro) \Y, 



2ro + 2i? 3MR 2iTmu'f 



ro 



+ 



+ ^{2rom^ + ^ro + Rm^){u^f 
2ro-4i? 3MR 2iTmu't' 



e(ro - r) 



'0 



'0 



2 27 9 ' 

+ -{-— R + 2rom'' + Rm' + -rQ){u'l'f 



Q{r-ro)\Y*^{eoAo), 



(E41) 



We note that only ff^'^ is discontinuous at r = ro. However, as mentioned before, the force depends only on Hq 
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and ii^S.RW .^^jiigji a,re continuous. We obtain the S-force as 
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Subtracting the above from the full force (E38), we find 
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We note that the Newtonian term, 2/i /rg, is precisely the correction to the force at 0{ijl ) when ro is the distance 
from the center of mass to the location of the particle. 
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If wc recall the fact that both //^""'^ and H^^^'"^ are discontinuous at r = rg and the gauge transformation from the 

Zerilli gauge to a RW gauge given by Eq. (E35) docs not change the discontinuity, while only '™ is discontinuous 
for the S part, we see that the RW gauge wc adopted to obtain the full force is different from the RW gauge for the 
S part obtained by the transformation from the harmonic gauge. Fortunately, however, because the force depends 
only on Hq (and K) for circular orbits, and its discontinuity structure at r = ro happens to be the same in both 
gauges, the resulting force (E43) turns out to contain no discontinuity. Furthermore, as discussed above, the correct 
Newtonian force is recovered at It is not clear if this desirable property holds because the orbit is circular or 

because only the IPN order correction is considered. If this happens to be no longer the case for general orbits, it 
will be necessary to find a gauge transformation that remedies the discrepancy. In any case, except for the correction 
at Newtonian order, the gauge ambiguity remains in the final result, and its resolution is left for future work. 



APPENDIX F: m-SUMMATION OF TENSOR HARMONICS 

In this Appendix, we summarize the formulas for summing over m-modes of the tensor harmonics for arbitrary £. 
Specifically, the m-sum we need to evaluate takes the form, 

e 

^ m^|y,„(7r/2,0)^ (Fl) 

where A'' is a non-negative integer. To perform the summation, we introduce the generating function, 

e 

Te{z)^ ^ e™1F,„(7r/2,0)|2. (F2) 
Then the sum (Fl) may be evaluated as lim d^Teiz). The above function is calculated as 

z— »0 

^dz) = ^e^F, (^i, 1; 1 - e'^^) , (F3) 

where 2F1 is the hypergeometric function. This can be easily expanded to an arbitrary order of z. For example, to 
0{z^), we have 

W = +[ 8 J 4!" 

£{£+ + - 51^ - 10^ + 8) ^ 1 ^6 ^ 



+ ^ ) 6! + I • 

In the cases of the vector and tensor harmonics, it is necessary to evaluate the m-sum of the form, 

e 

J2 m''\deYem{n/2,0)f . (F5) 
m=-e 

We introduce the generating function, 

e 

Ae{z)= J2 e'"19,F,„(7r/2,0)|2. (F6) 

m=-e 

This is expressed in terms of a hypergeometric function as 

2£+l^(,_,,, r(^+l/2)r(3/2) , f3 



A. W = ^e^'-'^^ - ^" .F, ^_ + 1; + e-J . (F7) 

The sum (F5) is evaluated by taking the derivatives of the above generating function. Expanding in powers of z, the 

m-sum (F5) is calculated as 



34 



+ 



+ 



16 



4! 



£{£ +1) {£-!){£ + 2) {5£'^ + 10£^-45£^-50£+m)\ 1 
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